Unravelling factors influencing demand for modern contraception and evaluating coverage progress since 2015 in Ethiopia, Kenya, and Nigeria: insights from multilevel and geostatistical modelling

Introduction The United Nations established the Sustainable Development Goals (SDGs) in 2015 to enhance global development. In this study, we examine an SDG indicator: the percentage of women aged 15–49 whose family planning needs are met by modern contraception (mDFPS). We evaluate both the factors influencing its coverage and its progress since 2015. Methods We used nationally representative surveys data (Demographic and Health Surveys (DHS) and Performance Monitoring for Action (PMA)) from Ethiopia, Kenya, and Nigeria. We assessed predictors of mDFPS. We also computed mDFPS coverage across countries and subnational areas, assessing coverage changes from the SDGs onset to the most recent period, using a Bayesian model-based geostatistical approach. We assessed whether the subnational areas exceeded the minimum recommended WHO mDFPS coverage of 75%. Results Varied individual and community-level determinants emerged, highlighting the countries’ uniqueness. Factors such as being part of a female-headed household, and low household wealth, lowered the odds of mDFPS, while rural-residence had low odds only in Ethiopia and Nigeria. The results indicate mDFPS stagnation in most administrative areas across the three countries. Geographic disparities persisted over time, favouring affluent regions. The predicted posterior proportion of mDFPS and exceedance probability (EP) for WHO target for Ethiopia was 39.85% (95% CI: [4.51, 83.01], EP = 0.08) in 2016 and 46.28% (95% CI: [7.15, 85.99], EP = 0.13) in 2019. In Kenya, the adjusted predicted proportion for 2014 was 30.19% (95% CI: [2.59, 80.24], EP = 0.06) and 44.16% (95%CI: [9.35, 80.24], EP = 0.13) in 2022. In Nigeria, the predicted posterior proportion of mDFPS was 17.91% (95% CI: [1.24, 61.29], EP = 0.00) in 2013, and it was 23.08% (95% CI: [1.80, 56.24], EP = 0.00) in 2018. None of the sub-national areas in Ethiopia and Nigeria exceeded the WHO target. While 9 out of 47 counties in Kenya in 2022 exceeded the WHO mDFPS target. Conclusion The study unveils demographic, geographic, and socioeconomic mDFPS disparities, signalling progress and stagnation across administrative areas. The findings offer policymakers and governments insights into targeting interventions for enhanced mDFPS coverage. Context-specific strategies can address local needs, aiding SDG attainment. Supplementary Information The online version contains supplementary material available at 10.1186/s12889-024-19387-9.


Introduction
The United Nations member states adopted the Agenda 2030 for Sustainable Development Goals (SDGs) in 2015.The agenda has 17 goals building on the Millennium Development Goals (MDGs) to stimulate action toward shared sustainable prosperity.The use of modern contraceptive methods falls under SDG 3.7, which states that "by 2030, countries should ensure universal access to sexual and reproductive health (SRH) care services for family planning, information, and education, and the integration of reproductive health into the national strategies and programmes" [47].Several indicators track this goal, one of which is indicator 3.7.1,which measures the percentage of women of reproductive age (15-49 years) who have their demand for family planning satisfied with modern contraceptive methods (mDFPS) [47].
Globally, an estimated 1.1 billion women of reproductive age needed family planning to postpone or avoid getting pregnant in 2020 [48].However, only 851 million had access to and used modern contraceptives, and an additional 85 million used traditional contraception methods.Modern contraception methods are defined as any product or medical procedure that prevents pregnancy from occurring as a result of sexual intercourse [25].Among those women who were using contraception in 2020, 90% of them were using a modern method [48].Modern contraceptive methods are considered more effective at preventing pregnancy than traditional methods, and many health facilities encourage their use (Sully et al. 2020).Sub-Saharan Africa (SSA) has the lowest proportion of modern contraceptive use among all women of childbearing age, with just 22.0% (95% Confidence interval: [21.8, 22.2]) [11].In 2020, among women who wanted to postpone pregnancy, only 55% were using modern contraceptives in SSA.In addition, most countries with lower than 50% mDFPS were from the region [14,48].
Women who can make informed decisions about their sexual relations, contraceptive use and reproductive health are more likely to use modern contraceptives to meet their family planning needs [6].In turn, meeting the demand for family planning with modern contraception methods empowers women to prevent unintended and high-risk pregnancies, thereby lowering the risk of maternal and under-five mortality [16,45].Additionally, modern contraception reduces the risk of early childbearing and allows women to achieve their educational goals.Thus, access to and use of modern contraception can improve women's health and educational opportunities, hence improving women's chances of paid labour employment [9,40].
Since the inception of the SDGs in 2015, it has been imperative that countries track their progress within different demographic, social, economic, and geographical groups to ensure that no one is left behind [5,47].Countries can identify disparities in access to modern contraceptives and other services by monitoring this progress.They can then use this information to address these disparities and ensure that all women of reproductive age can access the necessary information and services to achieve their desired family planning needs with modern contraception methods [47].
The current analysis focuses on Ethiopia, Kenya, and Nigeria, with Nigeria and Ethiopia being the most populous nations in Africa and Kenya ranking seventh [51].Kenya ranks among the top countries in SSA regarding the modern contraceptive use prevalence rate (mCPR) among married women, at 56.9%.Ethiopia, on the other hand, has a moderate mCPR rate of 40.5%, whereas Nigeria has one of the lowest prevalence rates among married women in SSA, at 12% [11,19,37].In terms of the proportion of mDFPS among married women, coverage was 64.9% in Ethiopia, 74.7% in Kenya, and 33.9% in Nigeria [19,37,48].
In addition, there is an unequal distribution of modern contraceptive use among various categories within these countries.Adolescents and young women (ages 15 to 24), economically disadvantaged women, women with limited education (below secondary school level), and rural residents have lower utilisation rate [2,10,12].Although the predictors of modern contraceptive utilisation have been investigated, there needs to be more literature regarding the predictors of mDFPS [5].Hypothesising that similar factors influence the adoption of modern contraception and mDFPS, this study seeks to investigate these factors and evaluate and compare the progress in mDFPS coverage since 2015 in the three countries.

Data sources
The study used data from Demographic and Health Surveys (DHS) and Performance Monitoring for Action (PMA) from three countries, namely Ethiopia, Nigeria and Kenya [19,30,31,36,54].The DHS is a nationally representative, cross-sectional survey that provides data on various health indicators, including family planning.PMA surveys cover a subset of the variables included in the DHS and are done in selected administrative areas in each country.DHS data, collected closer to the inception of the SDGs in 2015, was utilised.For more recent data where DHS was unavailable, the PMA survey data was utilised.Data from Ethiopia was obtained from the 2016 DHS survey and the 2019 PMA survey, excluding the 2019 mini-DHS survey due to the absence of necessary questions for calculating mDFPS [18,54].Kenya data was analysed using the DHS for 2014 and 2022 [30,31].Nigeria's data analysis relied on the 2013 and 2018 DHS surveys [36,37].
The outcome of interest was defined as a binary variable (Yes/No): where women 15-49 years old who needed family planning and were using modern contraception methods to satisfy that need were considered to have demand for family planning satisfied with modern contraception methods (Yes); the need for family planning was defined as fecund women who either needed to space or limit childbearing [5,20].By choosing to use mDFPS instead of just modern contraception prevalence, we were able to restrict the analysis to individuals with a demand for contraception, making it easier to track the progress of SDG indicator 3.7.1 [5,20].Modern contraceptive methods include pills, condoms (male and female), injectables, hormonal implants, patches, diaphragms, spermicidal agents (foam/gel) and emergency contraception.

Statistical analysis
We adjusted exploratory statistics for complex survey design [34].

Multilevel modelling
Using the latest standard DHS survey datasets for each country, we used a multilevel Bayesian logistic model to fit the data to explore factors associated with the outcome of interest.Multilevel models were used because of the hierarchical nature of the data's sampling framework (Eq. 1) [44].All Women from sampled households in enumeration areas (clusters) are selected for inclusion.The enumeration areas are, in turn, located in administrative areas (counties, regions or states).The independent variables included in the adjusted models were selected based on the literature, to allow within and between country comparison [20,44].The following variables were included in the multilevel model: place of residence (rural versus urban), age category, respondent's level of education, household head gender, wealth quintile, marital status and religion and parity.
Where the indices are defined for region i , cluster j and individual k .Index ji denotes that cluster j is nested within a region i.Similarly, index k: (j: i) denotes that individual k is nested within a cluster (whereupon the cluster is also nested within a region i ).The vector X denotes the covariates, and β sare the coefficients for the covariates.

Geostatistical modelling
A model-based Bayesian geostatistical logistic model to derive and assess the geographical variation in coverage for mDFPS (Eq.2) [32] was also used in the study.The variables used to adjust for confounding in these models were informed by the variables identified in the multilevel modelling stage and the literature.We included the variables in the analysis based on the availability of the spatial raster data of those variables for the countries of interest.The raster data included in the geostatistical models were women's education years, population density, probability of seeking care at health facilities, poverty rate and mean parity [52].
Where the vector X denotes covariates and β s are coefficients for the covariates conditional on the true prevalence P(x i ) at location S(x i ), i = 1, …, n, where the number of positive results is x i out of N i with a binomial distribution.S(.) is a spatial random effect that follows a zero-mean Gaussian process with the Matérn covariance function.
The geostatistical models were fit to two time points, at the start, closest to 2015 and end, the most recent, to allow comparison in the coverage of mDFPS between the two-time points to assess progress.Posterior samples were drawn based on the models to calculate sub-national-level mDFPS estimates and 95% credible intervals (CI) for each sub-national area.In addition, the predicted posterior samples were used to estimate the mean posterior change (MPC) of the outcome of interest between the years on which the geostatistical models were based to assess the change (progress) in coverage of mDFPS.Given the available data, we quantified how likely it was for the coverage of mDFPS to be above the recommended WHO minimum coverage of 75% [47] for each sub-national area (region, state or county), set as threshold t.This was achieved by quantification of the exceedance probability (EP) presented in Eq. 3.An EP close to 1 indicates that the prevalence of the outcome (1) was above t, EP close to 0 indicates that the prevalence was likely below t.An EP value of around 0.5 indicated that the prevalence was equally likely to be below or above t, hence the high uncertainty [32].Equation 3 was also used to calculate EPs for MPC between the two time periods to assist in assessing if there was an improvement in coverage between the two times.In this case, the threshold t in Eq. 3 was set to zero [8].
For each country, we generated coverage of the outcome, and EP maps at 5 by 5 km pixel-, sub-national-, and national levels to influence decision-making at various administrative levels.These maps show areas where coverage is the lowest or highest.These maps and tables would enable program implementers to appreciate where progress is lagging and to prioritise health interventions in areas that need them most while maintaining support where coverage is already good [29].R programming software version 4.2.1 was utilised for data management, and the Bayesian models were fitted using the INLA R package [32].

Results
The Ethiopian DHS 2016 survey had 5,312 women aged 15-49 who expressed a desire for family planning.The sample was drawn from 606 clusters.Similarly, the Ethiopian PMA 2019 survey included 3,396 women, sampled from 265 clusters.The 2014 Kenyan DHS had 7,840 women with a demand for family planning, selected from 1,549 clusters.The 2022 Kenyan DHS survey comprised 8,911 women with a demand for family planning, drawn from 1,662 clusters.The 2013 Nigerian DHS encompassed 11,464 women with data on the demand for family planning, sampled from 887 clusters, while the 2018 survey included 12,243 women with a demand for family planning, drawn from 1,376 clusters.

Respondent characteristics and mDFPS coverage
According to the latest standard DHS surveys of Ethiopia 2016, Kenya 2022 and Nigeria 2018, the majority of respondents were rural residents in Ethiopia (80.99%) and Kenya (59.75%) while in Nigeria rural residents were (46.27%).About 20% of respondents were aged 15-24, Ethiopia (23.31%),Kenya (24.12), and Nigeria (20.42%).Refer to Table 1 for more details on the characteristics of the respondents.

Multilevel modelling
The results of the multilevel Bayesian logistic model varied between countries.This analysis was based on the latest standard DHS from each of the three countries for easier comparison between countries.In Ethiopia and Nigeria, women residing in rural areas had significantly lower odds of mDFPS compared to women residing in urban areas.In Ethiopia, the OR was 0.59 (95% CI: [0.42, 0.82]), implying 41% lower odds of mDFPS among rural women.In Nigeria, the OR was 0.88 (95% CI: [0.78, 0.99]), signifying 12% decreased odds of mDFPS among rural women.But in Kenya, the OR of mDFPS was not significantly different between rural and urban residents; (OR: 0.91, 95% CI: [0.78, 1.06]).Detailed information is presented in Table 3.
In Ethiopia, women with secondary school or higher level of education did not have significantly different odds of mDFPS when compared to women with primary school education and no education (OR: 1.09, 95% CI: [0.85, 1.38]) and (OR: 1.14, 95% CI: [0.92, 1.41]).However, in Kenya, only women with no education had lower odds of mDFPS (OR: 0.40, 95% CI: [0.31, 0.50]) versus women with a secondary school or higher-level education.Similarly, in Nigeria, only women who did not have any education had lower odds of mDFPS (OR: 0.51, 95% CI: [0.44, 0.59]) versus women with a secondary of higher-level education (Table 3).These findings suggest that the relationship between education level and mDFPS varies across the three countries, with education being positively associated with higher mDFPS in Kenya and Nigeria but not in Ethiopia.
In terms of marital status, women who had never been in a union had higher odds of mDFPS compared to those who were in a union in Kenya (OR: 1.41, 95% CI: [1.16, 1.72]) and Nigeria (OR: 1.51, 95% CI: [1.21, 1.87]), but not in Ethiopia (odds ratio: 0.80, 95% CI: 0.53, 1.20).Furthermore, women who were no longer in a relationship had higher odds of mDFPS compared to those who were in a union: Ethiopia (OR: 3.  3).
In comparison with Muslim women, those who identified themselves as Orthodox, Roman Catholic or Protestants had higher odds of mDFPS.Women with a parity of 3-4 did not have significantly different odds of mDFPS to those with a parity of none in Ethiopia and Nigeria, while in Kenya women with a parity of none had lower odds of mDFPS (OR: 0.38, 95% CI: [0.29, 0.49]).Women with a parity of 1-2 had higher odds of mDFPS in Ethiopia (OR: 1.49, 95% CI: [1.22, 1.82]), while women with a parity of five or more had lower odds of mDFPS in Ethiopia (OR: 0.64, 95% CI: [0.53, 0.77]) and Kenya (OR: 0.78, 95% CI: [0.67, 0.90]).Furthermore, in all three countries, there were substantial residual variation at both the cluster level and the regional, county, or state level, after accounting for the confounding variables (Table 3).

Geostatistical modelling
The results of the Bayesian geostatistical logistic regression based on the latest survey are presented in Table 4.The implemented models were adjusted to account for  S3.
The predicted posterior samples of the proportion of mDFPS from the geostatistical models were used to map the adjusted mDFPS of the three countries between the two time periods (Figs. 1, 2 and 3).The adjusted predicted posterior proportion of mDFPS and EP for exceeding the 75% WHO target for Ethiopia was 39 The pattern of mDFPS coverage was similar between the two time periods for all three countries, with some exceptions (Figs. 1, 2 and 3).None of regions in Ethiopia had an EP of greater than or equal to 0.90, surpassing 75% WHO recommended coverage in either period.In Kenya, in 2014, seven counties (Embu, Kiambu, Kirinyaga, Machakos, Murang'a, Nairobi and Nyeri) had an EP of greater than or equal to 0.9 that they exceeded the 75% WHO recommended coverage, while in 2022 there were only seven (Embu, Kiambu, Kirinyaga, Machakos, Murang'a, Nyandarua, Nyeri), with an additional two counties (Nairobi and Tharaka-Nithi) that had an EP of above 0.8.In Nigeria, none of the states had an EP probability greater than or equal to 0.9 in either survey year in 2013 and 2018.
In addition, the geostatistical predictive posterior samples from the geostatistical models were used to calculate the mean posterior change (MPC) in the proportion of mDFPS between the two time periods for each included country.The MPC was calculated overall for the country and at the sub-national level (Table 5, Supplementary material Tables S1 & S2).An EP was also calculated, interpreted as the probability that the MPC between the first and the second period was greater than zero; in other words, to assess whether mDFPS coverage improved in the second period.The overall national level 95% CI of the absolute change in the proportion of mDFPS for the three countries were as follows: Ethiopia: 5.68% (95% CI: [-38.07,49.44], EP = 0.61); Kenya: 10.19% (95% CI: [-17.72,39.33], EP = 0.80) and Nigeria: 1.98% (95% CI: [-26.48, 34.49], EP = 0.58).
Kenya had the biggest absolute national change in the estimated percentage change of mDFPS coverage between the two periods.It also had a high probability that the absolute difference was greater than zero.The EP of the change of mDFPS being greater than zero was closer to 0.60 for Ethiopia and Nigeria.At the sub-national level, in Ethiopia, only one region, Gambela, had a high EP probability (0.90) that the change in mDFPS coverage proportion improved in the second period of 2019 compared to 2016 (Table 5).In Kenya, the counties that had an EP of 0.9 or more probability of improving in mDFPS were 22 counties (Baringo, Bomet, Bungoma, Elgeyo-Marakwet, Homa Bay, Isiolo, Kajiado, Kericho, Kilifi, Kisumu, Laikipia, Mandera, Migori, Nakuru, Narok, Nyandarua, Taita Taveta, Trans Nzoia, Turkana, Vihiga, Wajir, West Pokot) between the year 2014 compared to 2022 (Table S1).While in Nigeria, the states with an EP of 0.9 or more probability of improving in mDFPS were 9 (Adamawa, Buachi, Gombe, Jigawa, Kano, Katsina, Kebbi, Sokoto and Zamfara) between the years 2013 compared to 2018 (Table S2).Conversely, the areas that had a very low EP probability ( < = 0.1) of improvement in coverage of mDFPS were likely to have had a reversal in coverage of mDFPS in the later period (Table 5 and Supplementary material Tables S1 & S2).

Discussion
The results from multilevel modelling varied among the three countries, except for a few determinants, emphasising the uniqueness of each country.For instance, in Ethiopia, women without an education or with primary-level education had similar odds of mDFPS as those with secondary or higher-level education.In Kenya and Nigeria, women without an education had lower odds of mDFPS.This suggests that education may not contribute to mDFPS uptake inequality in Ethiopia compared to Kenya and Nigeria.In all three countries, women from well-todo households had higher odds of mDFPS, highlighting that household wealth remains an important determinant of inequality [1,4].Muslims consistently had lower odds of mDFPS than women from other denominations; this provides an opportunity for tailored interventions to improve this population's access and use of modern contraception [6].
In all three countries, geographic variations in the proportion of mDFPS were observed.The pattern of geographic variation was similar for the first and second periods, with a higher proportion of mDFPS seen in the relatively affluent sub-national areas.The analysis results investigating whether there was an improvement in the proportion of mDFPS, show that both Ethiopia and Nigeria had modest improvements.In contrast, Kenya increased its coverage of mDFPS between the two periods; 22 out of the 47 states improved coverage of mDFPS over the two periods, 2014 to 2022, with 9 counties likely to have exceeded the minimum WHO recommendation threshold in the second period, 2022.In the sub-national analysis in Ethiopia, Gambela region (out of 12 regions) was identified as the only region experienced a significant increase in its mDFPS coverage, over the two periods, 2016 to 2019.In Nigeria, 9 (out of 36 states) were identified as having had an increase in mDFPS coverage from 2013 to 2018, all from the northern region.In Ethiopia and Nigeria, most administrative areas had only slight improvements or stagnated in mDFPS coverage.It should be noted that even in regions where improvements were registered, the level of coverage of mDFPS was still below the WHO target.This can be seen on the exceedance maps of the WHO recommended 75% coverage that barely changed between the two periods.
The current analysis results indicate that demographic, geographic, and socioeconomic factors are important in determining who has greater coverage of mDFPS.These findings are consistent with previous studies that have analysed data on mDFPS [5,44].The geospatial analysis identified areas that still need to meet the WHO mDFPS coverage target.The geospatial analysis of changes in the coverage of mDFPS between the two time periods identified regions that have experienced either an increase or a decline and those that have been stagnant.These findings can guide targeted interventions that can address the unique needs of local communities instead of generalised interventions that local evidence does not drive [3,20,39,53].
Administrative areas with declining or stagnant mDFPS coverage should be prioritised for interventions to improve the situation.In areas that have shown progress, efforts should be sustained to maintain and accelerate the gains made [45].Even though we are advocating for interventions to increase the level of coverage of mDFPS, the rights of the individual and the couples and communities must be respected as enshrined in the 1994 International Conference on Population and Development (ICPD) conference [22].It is important to acknowledge that some women may not start using modern contraceptive methods soon, even if they have an unmet need for modern contraception methods due to reasons such as perceived low risk of pregnancy or perceived cultural, social, or health concerns and preference for traditional methods [15,42].On the other hand, it's worth noting that a portion of the unmet need for contraceptives may reflect a desire for modern contraception use, making these individuals more likely to adopt it when made available and accessible [15,17].
Access and utilisation of modern contraceptive products can also be affected by conflicts and climatic emergencies [35,46].Women living in areas with ongoing conflict or climatic catastrophic events have reduced usage of modern contraception [35,46].In general, conflicts and climate change events can pose significant challenges to countries in meeting the demand for family planning services among women.When such events occur, the health system may be overwhelmed and unable to provide the necessary services that enable women to access and use modern contraception methods [46].These challenges can have a particularly profound negative effect in regions and countries with already vulnerable healthcare infrastructure, like most countries in the SSA region [21,48].
In Ethiopia, the Tigray region (and other regions such as Amhara, Afar and Oromia) has been at war since November 2020 until the time of writing this paper.Since the data analysed in this work was from a pre-conflict period, the situation is likely to have changed [21].The Afar, Oromia, Somali and Harari regions are prone to drought and were also identified as having low mDFPS coverage in this analysis [7,24,28].Similarly, in Kenya, most counties (Garisa, Isiolo, Mandera, Marsabit, Tana River, Samburu, Turkana and Wajir) with a history of drought were also identified as having a low coverage of mDFPS.When communities live in resource-constrained conditions, they prioritise survival and hence have low utilisation of health services [33,49].
Furthermore, the northern Kenyan counties also face conflicts as communities are forced to fend for scarce resources.This results in displaced communities, further exacerbating the issues of health service utilisation [28].To promote the use of modern contraception methods in such areas, it is necessary to take a holistic approach that addresses the challenges caused by conflict and climate change shocks [41].
In Nigeria, the northern states have been at war since 2009 [38]; our analysis shows that most of the northern states have low coverage of mDFPS.The region is also predominantly Muslim [43].Our analysis and work from elsewhere show that Muslim women have lower odds of modern contraception use than women from Christianbased religions in the SSA [6].In addition, misconceptions about modern contraception methods discourage women from using these methods in northern Nigeria.Some documented misconceptions include the belief that modern contraception can affect a woman's ability to bear children and can cause cancer [26].Work being done by humanitarian organisations can be alluded to as having helped to improve the level of mDFPS between the two time periods in the northern states, but there is more work that needs to be done to improve healthcare services in the region [27,50].Overall, Nigeria is one of the countries in SSA with a low modern contraception prevalence rate.In 2017, the Healthy Policy Plus initiative identified insufficient domestic funding as one of the challenges affecting Nigeria's family planning programmes [23].

Study strengths and limitations
The study had several strengths.Firstly, we employed a multilevel modelling approach to identify the disparities in mDFPS coverage driven by demographic, geographic, and socio-economic factors.This modelling process helped us identify the cluster and regional-level residual variance, indicating the presence of spatial variation.To account for this, we performed a spatial analysis using geostatistical models.This enabled us to calculate the adjusted proportion of mDFPS for each administrative area and determine the exceedance probability of each area, highlighting areas which exceeded the WHO mDFPS coverage minimum target of 75%.Additionally, our analysis of two time periods, which corresponded to a before and after period, allowed us to assess the progress made and identify areas that have stalled or retrogressed in their mDFPS coverage.The methods used in this study can be replicated for other SDG indicators and other important indicators to identify underserved communities, assisting policymakers and governments in directing their efforts more effectively.Our research also had certain limitations.Despite our efforts, we could not find nationally representative data from 2015, when the UN established the SDG goals.Instead, we used the data that was closest to 2015.Similarly, for the latest year, we relied on the most recent data available for these countries.Thus, the time between the assessment years was different among the included countries, affecting the comparability of the results between the three countries.A complete spatialtemporal approach was impossible due to limited available data, which only spanned two time points.It is hard from this analysis to tell the actual reasons why the proportion of mDFPS reduced or improved in the areas that it did; hence, further studies are required in these areas.The survey responses may have been affected by social desirability bias, in which respondents provided answers that they believed were more favourable.It is also possible that the survey participants might have had recall bias; where they could have not accurately remembered the information they were being asked.Also, methodologies used between the two survey periods could have had changed, making direct comparisons challenging.Caution is necessary when generalizing these results to other countries, as their contexts may differ from the countries included in this analysis.The predicted posterior 95% CI's were wide, but the uncertainty is comparable to DHS spatial analysis guidance recommendation [13].Despite these limitations, our analysis provides valuable evidence and fills an essential gap in the available literature.

Conclusion
Our analysis utilised multilevel modelling, geospatial analysis and comparison of data from three countries over two time periods.The results showed the presence of demographic, geographic and socio-economic disparities in the mDFPS, and that some areas have made progress while others have retrogressed but the majority have remained stagnant.The majority of sub-national level areas have mDFPS that are below the WHO recommendation of 75%.Our findings will be important in providing valuable insights for policymakers and governments on how they can effectively target their interventions to improve uptake of modern contraception among women who have a demand for it.By providing much-needed evidence, this analysis contributes to countries' efforts to assess their progress in meeting the SDG indicator 3.7.1.

Fig. 1
Fig. 1 Proportion of demand for family planning satisfied with modern contraception methods (mDFPS) for Ethiopia.Estimates based on adjusted Bayesian geostatistical models.Top left panel: Map of the predicted proportion of demand for family planning satisfied with modern contraception methods (mDFPS) in Ethiopia in 2016 (DHS).Top right panel: The likelihood or certainty that the estimated mDFPS exceeds the 75% threshold (SDG target of mDFPS).Bottom left panel: Map of the predicted proportion of mDFPS in Ethiopia in 2019 (PMA).Bottom right panel: The likelihood or certainty that the estimated mDFPS exceeds the 75% threshold (SDG target of mDFPS)

Fig. 2
Fig. 2 Proportion of demand for family planning satisfied with modern contraception methods (mDFPS) for Kenya.Estimates based on adjusted Bayesian geostatistical models.Top left panel: Map of the predicted proportion of demand for family planning satisfied with modern contraception methods (mDFPS) in Kenya in 2014 (DSHS).Top right panel: The likelihood or certainty that the estimated mDFPS exceeds the 75% threshold (SDG target of mDFPS) .Bottom left panel: Map of the predicted proportion of mDFPS in Kenya in 2022 (DHS).Bottom right panel: The likelihood or certainty that the estimated mDFPS exceeds the 75% threshold (SDG target of mDFPS).

Fig. 3
Fig. 3 Proportion of demand for family planning satisfied with modern contraception methods (mDFPS) for Nigeria.Estimates based on adjusted Bayesian geostatistical models.Top left panel: Map of the predicted proportion of demand for family planning satisfied with modern contraception methods (mDFPS) in Ethiopia in 2013 (DHS).Top right panel: The likelihood or certainty that the estimated mDFPS exceeds the 75% threshold (SDG target of mDFPS) .Bottom left panel: Map of the predicted proportion of mDFPS in Ethiopia in 2018 (DHS).Bottom right panel: The likelihood or certainty that the estimated mDFPS exceeds the 75% threshold (SDG target of mDFPS).

Table 2
Coverage of demand for family planning satisfied with modern methods (mDFPS) by characteristics of respondents based on latest DHS data.Both counts and percentages were weighted to adjust for survey design.Additional details are in tables S4 to S6

Table 3
Table of adjusted odds ratios (aOR) for demand for family planning satisfied by modern contraception methods (mDFPS) in Ethiopia, Kenya and Nigeria using bayesian logistic regression (eq. 1) based on latest DHS data

Table 4
Adjusted odds ratios (aOR) for demand for family planning satisfied with modern contraception methods in Ethiopia (PMA), Kenya (DHS) and Nigeria (DHS) based on the latest health surveillance survey using bayesian geostatistical logistic regression (eq.2) 1Percentage of people living below US $1-dollar variable was not available for Ethiopia, PMA Performance Monitoring for Action, DHS Demographic Health Survey

Table 5
Table showing Ethiopia's regional proportion of mDFPS (2016-2019) with 95% credible intervals (CI), mean posterior change (MPC), percentage change, and exceedance probability of increase of greater than zero of mDFPS from 2016 to 2019.Estimates based on adjusted bayesian geostatistical models